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Abstract 

We obtain the exact nontopological soliton lattice solutions of the Associated Lame equation in 
different parameter regimes and compute the corresponding energy for each of these solutions. We 
show that in specific limits these solutions give rise to nontopological (pulse-like) single solitons, as 
well as to different types of topological (kink-like) single soliton solutions of the Associated Lame 
equation. Following Manton, we also compute, as an illustration, the asymptotic interaction energy 
between these soliton solutions in one particular case. Finally, in specific limits, we deduce the 
soliton lattices, as well as the topological single soliton solutions of the Lame equation, and also 
the sine-Gordon soliton solution. 

PACS numbers: 03.50.-z, 05.45.Yv, ll.lO.Lm 
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I. INTRODUCTION 



Over the years, extensive research has been carried out seeking the exact sohton solutions 
of both periodic (e.g., sine-Gordon, double sine-Gordon) and nonperiodic (e.g., 0^,0^) field 
theory models. For example, the exactly solvable sine-Gordon (SG) equation |1| and its 
quasi-exactly solvable (QES) partner, i.e., the double sine-Gordon equation (DSG), have 
exact single soliton 2, 0] as well as soliton lattice 4] solutions. 

There have been some advances in the study of the hyperbolic analogues of these prob- 
lems; the exactly solvable hyperbolic analogue of the SG equation is the sine-hyperbolic 
Gordon (ShG) equation ^j. This potential has only one minimum and thus does not sup- 
port (topological) soliton solutions. The hyperbolic analogue of the DSG equation is the 
double sine-hyperbolic Gordon (DShG) equation^which is a QES double-well potential with 
exact single soliton and soliton lattice solutions p. 

However, not much is known regarding the elliptic analogues of these problems. The 
elliptic 'generalization' of the SG is the Lame equation Q but, as far as we are aware 
of, its single soliton and soliton lattice solutions have not been worked out yet. One would 
surmise that the elliptic analogue of the DSG equation is the Associated Lame (AL) equation 
0, y, Q. Nevertheless, we find below that this is not the case. 

Generally speaking, the Schrodinger equations with periodic potentials belong to a class 

energy spectrum of 

the system. For both the SG and DShG potentials one could calculate the exact statistical 
and thermodynamical properties 0, through a knowledge of the spectral band structure 

utions, respectively. Both AL and Lame are also periodic 
y, [l^ . Although our goal in the present paper is only to 
obtain their exact single soliton and soliton lattice solutions, this will serve as a background 
for future statistical mechanics studies of the AL and Lame elliptic potentials. In particular, 
we obtain here the soliton lattice solutions of the AL equation in different parameter regimes 
corresponding to different shapes of the AL potential. The advantage of this approach is that 
the corresponding single pulse and topological soliton solutions of the AL equation can be 
immediately obtained from the lattice solution in suitable limits. Furthermore, the soliton 
lattices and topological single soliton solutions of the Lame equation are easily obtained by 
taking a different appropriate limit. Besides, the asymptotic interaction between the solitons 



and their solitonic and phononic so^ 
potentials with a band structure [? 
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121. 



can also be readily obtained by using Manton's formalism 

Our results can be summarized as follows: We show that there are six different soliton 
lattice solutions of the AL equation for different values of the parameters. From these 
six solutions we obtain the corresponding AL single soliton solutions. It turns out that 
while in five cases these are topological (kink-like) single solitons, in one case we have a 
nontopological (pulse-like) single soliton solution. From the AL soliton lattice solutions, by 
taking appropriate limits, we also obtain the corresponding Lame soliton lattices and single 
soliton solutions. 

There are many physical contexts in which the Lame equation arises, such as bond-order 
and charge density wave systems [13], nonlinear elasticity 14], and other contexts, e.g., 
phase slips in superconductors and magnetoelastic interaction on curved surfaces. The AL 
equation could arise in the above systems in the presence of external - electric, magnetic or 
stress - fields. 

The plan of the paper is as follows. In Sec. II we discuss some salient features of the 
AL potential. In Sec. Ill we obtain the nontopological soliton lattice solutions of the AL 
equation in different regimes in the parameter space and also compute the corresponding 
energies. In Sec. IV we obtain the topological and pulse single soliton solutions of the AL 
equation by taking the appropriate limits of the various nontopological solutions obtained 
in the previous section. It is worth emphasizing here that in one special case, we obtain 
two different kinds of topological solutions of the AL equation. As an illustration, in one 
particular case, we also estimate the asymptotic interaction between these single solitons. 
(In all the other cases one can follow exactly the same procedure, and thus we will not 
consider these here.) As a crosscheck on our results, we recover in appropriate limit, the 
sine-Gordon soliton solution. In Sec. V we show that by taking an appropriate limit, we can 
also obtain the kink lattice and the topological single soliton solutions of the Lame equation. 
Finally, in Sec. VI we summarize the results obtained in this paper and indicate some open 
problems. 
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II. THE ASSOCIATED LAME POTENTIAL 



Consider the following family of periodic potentials labeled by a pair of real parameters 
(p, q): 

K4l(0, k) = pPsn^i^, k) + qk^^^^^^^ + C = p^sn^icl), k) + sn^cj) + K{k), k) + C 

dn (0, k) 

(1) 

that are called Associated Lame potentials (since the corresponding Schrodinger equation is 
called the Associated Lame equation) 0,0,13]. Here sn(0, k) and cn(0, k) are respectively 
the sine and cosine amplitude Jacobi elliptic functions of real modulus (0 ^ k ^1) and 
period AKiJi); dn(0, k) is the 5-amplitude Jacobi elliptic function of modulus k and period 
2 K{k); and K{k) denotes the complete elliptic integral of the first kind, see jj^lffij]. We will 
choose the constant C in the potential so that the absolute minimum of the potential-with 
respect to 0-equals zero, Vmin = 0. The potential is periodic with period 2K{k), except in 
the limit p = q when the period reduces to K(k), as it is clear from Eq. (^. 

The case q = corresponds to the standard Lame potential Vl{4>, k). There are two 
more cases in which the potential reduces to the Lame potential. Namely, when p = 
and also when p = g, as shown below. Therefore in all the calculations below we will 
admit that p ^ and p ^ q] the results for the standard Lame potential will be recovered 
in the limit q ^0, and we refer to these results whenever the potential reduces to the 
standard Lame one. From a physical point of view, if one thinks of a Lame potential (p, 0) 
as due to a one-dimensional array of atoms with spacing 2K{k) and 'strength' p, then the 
Associated Lame potential {p, q) results from two alternating types of atoms spaced by K{k) 
with 'strengths' p and g, respectively. If the two types of atoms are identical (which makes 
p = q), one expects a potential of period K{k). 

In addition, when /c — 0, and |g| oo, so that \p\ k"^ ^ P = finite and \q\ k"^ ^ Q 
= finite, the potential ((1} reduces to the sine-Gordon potential, 

VsG = {P sign p - Q sign q) sin^ , (2) 

and the results we present below reduce to the well-known ones for the sine-Gordon potential. 

It is also worth noting that under the transformation 0^0 + K{k) the AL potential 
Val(p, (?) goes over into Val{<1,p) and hence for p > 0,g > 0, as well as for p < 0,q < 0, 
without loss of generality, we shall always consider the case of > g^. Further, instead of 



considering both the possibihties of p > 0, g < and j9 < 0, g > 0, it suffices to consider just 
the case of p > 0, g < 0, but now can be bigger as well as smaller than q\ 

Finally, consider the case of p = g, when, by using the Landen transform jl5| and choosing 
the constant C = —pk^, the AL potential ([[J can be written as 



l + k' 



(3) 



where k' = ^/l — k^ is the complementary elliptic modulus. For the rescaled field = 
(1 + k')(j) and in the rescaled space coordinate a; = (1 + k')x the field equations [see next 
section, Eq. ((H))] will remain the same as that for the simple Lame potential. Note that the 
p = q case cannot be obtained as a limit of g p (g 7^ p), since-as already mentioned-when 
p = q the periodicity of the potential (|T)) is K{k), while for p 7^ g it is 2K{k), see also j^. 



III. SOLITON LATTICE SOLUTIONS OF THE AL POTENTIAL 



For the scalar field ( 
differential equation 

In the stationary case 



(x, t) the dynamics is described by the second-order hyperbolic 



(x), this reduces simply to 

^ ^ BVal 
dx"^ d(h ' 



(4) 



(5) 



that can be easily integrated, at least formally, by quadratures. The time-dependent solu- 
tions are immediately obtained from here by Lorentz boosting. The physically meaningful 
solutions (bounded at x — > ± 00) are given by 



± V2\ 



X 



Xo) 



VaS, k) - 



(6) 



where xq and A = \p\ k a are the two (suitably chosen) integration constants. Of course, 

the absolute minimum/maximum of the po- 
tential energy with respect to (f). In fact, in view of the condition Val{4>, k) — A^ ^ 0, the 
actual value of determines the appropriate integration domain in Eq. (jH)). 

Depending on the values of the parameters p and g, the potential Val(0, k) has different 
behaviors in one period ^ < 2K{k), and hence the nature of the solutions in Eq. (jH)) 



is also different. There are three cases to be considered separately, namely (I) when p > 
and g ^ 0, (II) when p < and g ^ 0, and, finally, (III) when p > and g ^ 0. Note 
that in what follows we shall use the shorthand notation T = \/\q\/\p\. 



A. Case I : p > and g ^ 

As explained above, in this case it is sufficient to consider ^ F < 1. Note that the 
case r = corresponds to the standard Lame potential. Depending on the value of F, the 
potential Val{4>, k) can have different behaviors in one period ^ < 2K{k). 



Case /._?.• ^ r ^ k' 



Here ^ k' = \J\ — k'^ ^ 1 is the complementary elliptic modulus. The potential 
Val(0) has only one minimum, = (at = 0), and one maximum, Vmax = pk^ {1 — 

r^) [at (f) = K{k)], in one period [note that we have chosen C = — pfc^F^ in Eq. (Q) in 
order to have Vmin = 0]. The plot of the potential Val as a function of is given in Fig. 1(a), 
where the solid and the dashed horizontal lines correspond to two choices of the parameter 
A^. Under the change of variable z = sn^(0, k), Eq. (jH)) can be rewritten as 

(7) 



2^j2pk^ 
0, k), Xq 



X 



± 



Iz v/(zi - z){l - z) {z - Z2) {z - 0) 

where z = sn^(0, k), Xq = for 0o = K{k) (by choice) and = pk'^a? with ^ < 
[1 — F^). The value of a? determines the limits of the integration interval in the above 



equation, through the condition V(0i,2) 
given by: 



pk a = 0; one finds that zi^2 



sn 



>i,2j are 



1 - V'k"' + k^ 
2¥ 



^1,2 



1 ± 



4P 



[1 - T^k'^ + a2A;2)2 



(8) 



and zi>l>z^Z2>0. Then the integral in (|7j) can be evaluated using the formula 
3.147(5) of [l^ and one is finally led to the nontopological soliton lattice solution 

[{zi - Z2) - zi{l - Z2) sn2(y, t)\ 



sn (0, fc) 



[(zi - Z2) - (1 - ^2)sn2(y, i)] 



(9) 
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Here y = a/2 p A;"' (2:1 — Z2) x and the period of the lattice 2 L = 2 K{t)/ \j2pk'^ {zi — Z2) 
is controlled by the modulus 

^^ 1/2 



< t 



zi(l 



Z2 



^ 1. 



(10) 



Note that 



Zl - Z2 

oscillates between K{k) and sn~^(yz2) [respectively, — K{k) and 



— sn~^(yi2)], i.e., it is indeed a nontopological solution [it does not connect two adja- 
cent degenerate minima of the potential, that are separated by a distance-in the (j) space-of 
2K{ky\. In Fig. 1(b) we give a plot of the solution (i.e., the field as a function of x) with 
solid and dashed curves corresponding to the two choices of as in Fig. 1(a), while the 
dotted curve represents the topological single soliton solution (corresponding to = 0)-see 
Sec. IV. We also note that there are no "cusps" in the actual profiles of the fields; indeed, 
in Fig. 1(b) and in the subsequent figures we have plotted modulo 2K{k) the absolute value 
of the actual field 0(x). Specifically, for solutions that touch K{k) such as in Fig. 1(b), the 
actual field is (f){x) in [0, 2L], 2K{k) — 0(x) in [2L, 4L], 0(x) in [4L, 6L], and so on. Similarly, 
for solutions that touch field value such as in Fig. 3(b), the actual field is (f){x) in [0, 2L], 
— 0(x) in [2L,4L], 0(x) in [4L, 6L], and so on. 

One can compute the energy corresponding to a period 2 L of this soliton lattice: 



SL 



L 



dx"^ 



+ VaM 



4 / VAMdx -2LA^ = Apk'^ [ 

Jo Jo 

n^[0(a;), k] 



sn 



L 



10 ldn'[(j){x), k] 
After some algebraic manipulations, using the solution 



dx = j ^ [21^al(0) - A^] dx 
'^^■^'x)^ k] dx 
l \ dx - 2Lp2A;2a^ 
one obtains: 



(11) 



E 



SL 



2p 



Zl - Z2 



Zl k 



'2 



ZiKit) 



Kit) 



1^1 



i)n 



Z2 



Zl 



Z2 



1 



n 



;i - Z2) (1 - Zl k' 



1 — Zik"^ 1 — Zik'^ \ {zi — Z2) k"^ 

where n(zo,i^) denotes the complete elliptic integral of the third kind [151. II 



(12) 
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FIG. 1: Case I.l : ^ F ^ k' . (a) The shape of the potential Val{4>)- The sohd and the dashed 
horizontal lines correspond to two choices of the integration parameter A?, (b) The soliton lattice 
solution [i.e., the field (t>{x)] with solid and dashed curves corresponding to the two values of in 
(a). The x-axis is not labeled because the two lattices have different spatial periods. The dotted 
curve represents the topological single soliton solution corresponding to = (for the sake of 
clarity we have displaced the origin of the x-axis for this curve). 

Case 1.2 : k' < V < I 

The potential 1/^^(0) has now two minima: Vmin{<f> = 0) = and a local 
minimum Vj^^^[(f) — K{k)] — pk^{l — F^), and two symmetric maxima around 



= K{k), namely V^ax = pk'^ [(1 - rk')/k]'^ for = sn-^ ^(1 - rk')/k^ , 2K{k) 
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sn ^ 



v/(l - Tk')/k^ . The plot of the potential Val{(P) function of is given in 
Fig. 2(a), where the solid and the dashed horizontal lines correspond to the two choices of 
the integration constant as explained below. Note that in this case it makes no sense 
to consider either the limit of the Lame potential (i.e., F — 0), or that of the sine-Gordon 
potential (i.e., k ^ 0, p ^ P, q Q). There are two possible situations, depending 

on the value oi = p k^ a?. 

Case I.2(i) : Knm ^ =pk'^a^ < Vj^^^, i.e., ^ < 1 - F^. 
One recovers the same soliton lattice solution as above, Eq. Q, with the same modulus t, 
Eq. (fTUj) . and the same energy per period of the lattice, Eq. (fT^ . 

Case I.2(ii) : V^^,^ ^ A^=pea'^ < Vmax, i.e., 1 - F^ ^ < [(1 - Vk')/kf. 
Then the integral (jB)) becomes 

2^2pk^x = ± [ I ^ ^ ^ , (13) 



Z2 



V/(1 - Z) {Z, - ~Z)CZ - Z2) {z - 0) 



with z = sn^(0, k) and 1 > Zi ^ z > Z2 [-21,2 have the same expressions as above, Eq. (jHJ]. 
One can evaluate this integral using the formula 3.147(4) of [l5| and obtain the following 
nontopological soliton lattice solution 

^A<P, k) = ^ , (14) 

Zi — [Zi — Z2)sn^[w, n) 

with w = a/2p A;^ Zi (1 — Z2) x, and the modulus n given by 



therefore, the soliton lattice has a spatial period 2L = 2 K{n) / ^J2pk'^ zi{l — Z2). Note 
that varies between sn~'^{^/zi) and sn~-^(i/i2) [respectively, — sn~-^(i/ii) and — sn"^(i/52)] 
when X varies over one spatial period: indeed, this solution is nontopological. 

In Case I.2(ii) in order to compute the energy for one spatial period 2L, it is useful 
to shift the potential by — l^^j„ = —pk'^{l — F^), so that the minimum of the potential 
explored by the soliton lattice is equal to zero. One finds: 

EsL = ^ / VAMdx - 2LVi„ - 2Lp'ea' = 
J 



2p2 



Zl (1 - Z2) 



Zl ' / 1 — k'^ Z2 \Zi(l — k'^ Z2 



(16) 



The plot of the solutions in Case 1.2 is given in Fig. 2(b), where the solid line corresponds 
to the nontopological soliton lattice solution of Case I.2(i), while the dashed line corresponds 
to the nontopological soliton lattice solution of Case 1.2 (ii). Note that in the appropriate 
limits (see Sec. IV) these two types of solutions give birth to a topological [Case I.2(i)] and 
a nontopological single soliton [Case I.2(ii)] solution, respectively. 



B. Case II : p < and q ^ 

As explained above, here again it suffices to focus only on the domain ^ F < 1, and 
we must distinguish between two different behaviors of the potential, depending on the value 
of F. 



Case III .■ ^ r ^ /c' 



The potential Val{4', k) has only two extrema in [0, 2K{k)), namely an absolute maxi- 
mum Vmax = \p\ k"^ (1 — r^) for = 0, and an absolute minimum Vmin = for = K{k) 
[with the choice of the shift C = \p\ k"^ in Eq. (J])]. A plot of the potential is given in 
Fig. 3(a). 

One can repeat the integration scheme described in the preceding case. In particular, 
considering z = sn^(0, k), Eq. (jH)) becomes: 

dz 



2 v/2[prPa; = ± 



V(zi - 5) (1 - z) {Z2 - z) (5 - 0) ' 



where this time 



^1,2 



k^ - F2 



a2 e 



2F 



1 ± Wl 



4Pri - F2 - a^] 



(17) 



(18) 



(1 + F - F2 - a2 A;2)2 

Zi > 1 ^ Z9. > z > 0, and ^ < (1 — F^). This integral can be evaluated using formula 
3.147(2) of [15[ and one obtains the following nontopological soliton lattice: 

zi Z2sn^{y, r) 



sn (0, k) 



(19) 



{zi - Z2) + ^2sn2(?/, ^) • 
Here y = ^2 \p\ k^ ( Zi — Z2) X and the modulus r that controls the density of solitons in 
the lattice is given by 

0<r=,/^5LZil^l. (20) 



Zl - Z2 
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A' 



A' 



(a) 


1 








/ 1 \ 







K(k) 



2K(k) 



K(k) 




FIG. 2: Case 1.2 : A;' < T < 1. (a) The shape of the potential Val{<^). The sohd and the 
dashed horizontal lines correspond to two choices of the integration parameter A^, namely, Case 
I.2(i) : Vmin ^ = pk'^a'^ < F^„, i.e., ^ < 1 - (solid line); and Case I.2(ii) : 
^Ln ^ = Pk'^a^ < Vmax, 1-6., 1 - ^ < [(1 - rk')/k]'^ (dashed line), (b) The soliton 
lattice solution [i.e., the field (f){x)] with solid and dashed curves corresponding to the two values 
of A'^ in (a), i.e., respectively, to Case I.2(i) and Case I.2(ii). The x-axis is not labeled because the 
two lattices have different periodicities. 

The period of the lattice is 2L = 2 K (r) / ^/ 2 \p\ {zi — Z2) and one notices that (j) varies 
between and sn~^(y^) [respectively, and — sn~^(y^) ] in one period. A plot of the 
solution is given in Fig. 3(b). 
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K(k) 




2K(k) 2L 4L 

X 



6l9 



FIG. 3: Case II.l : ^ T ^ k'. (a) The shape of the potential Val(0). The solid horizontal line 
corresponds to the integration parameter A^. (b) The soliton lattice solution [i.e., the field 4>{x)] 
corresponding to the value of in (a) . 

One can compute the energy for one period 2 L of this lattice, 

Z2 



EsL = 2 





p 








zi { 



-(^1 - 1) K{r) + ;2i n 



Z\ — Zi 



1 - ^1 



1 - ^1 A;2 



z\ - Z2 



a^K(r) 



(21) 
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Case II.2 : k' < T < 1 



In this case, the potential Val{4>, k) has two maxima and two degenerate minima in 
one period ^ < 2K{k). Namely, an absolute maximum, V^ax = \p\ (1 — T k'Y (for 
= 0); a relative maximum V^^^ = \p\ (T — k'Y [for = K{k)\. and two absolute minima 



Vjnin = situated symmetrically around K{k) [for = sn ^a/(1 — F k')/k'^, 2K{k) 



?,'nr^^J{l — Vk')/^^]. Note that we used a shift C = (1 — 2 F A;' + F^) in the expression 
(H} of the potential. A plot of the potential as a function of is given in Fig. 4(a). 

The very existence of adjacent degenerate minima of the potential separated by different 
barriers to the right and to the left (on the axis) implies that in this case one will have 
two different soliton lattices (and, correspondingly, two different topological single soliton 
solutions). This is similar to the DSG case 

Baa. 

We follow the same type of integration procedure as that described above. 
Case II.2(i) : When Knm ^ = \p\ ^ < 1/^^^, i.e., ^ < (F - k'f/k\ Eq. ®, 
under the change of variable z = sn^(0, k), leads to two different types of solutions. 
Solution 1 : for < z ^ Z2 < Zi < 1, one has: 



z 



dz 



2 ^2 bl k* X = ± , , (22) 

Jo v/(l - Z) {Z, - Z) {Z2 - 5) (z - 0) ' 



where now: 



_ 2{l -Tk')-a^k^ ( A[{l-Vk'f -a'k^] . 

^^'^ - 2P I V [2{l - Tk') - a^k^f ^ ■ 



Equation (j22|l can be integrated using formula 3.147(2) of thus obtaining: 



sn' 0, k) = \ ^ ; ^ , 24 

(1 - Z2) + Z2S10.\W, S) 



with w = ^y2\p\k^ zi{l — Z2) X and the modulus 



It is a nontopological soliton lattice, with oscillating between and sn~^(y^) [respectively, 
and — sn~^(yi^)]. Its energy for one period 2L = 2 K{s)/ ^y2\p\ k^ Zi {1 — Z2) is given 
by: 



EsL 



2\p\ 



Z\ (1 - Z2) 



(F - k'Y a 
P 1 



Kis) + n(--^, . 

i — Z2 

- F^n(-^,.)^ (26) 
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Solution 2 : for < Z2 < Zl ^ z < 1 one obtains from Eq. (0): 

2 v/2H^ x = ± = , (27) 

Jz V(l - Z) (5 - ^i) (Z - Z2) (z - 0) 

with 2; = sn^0, k) and Zi^2 given by Eq. (PHj) . This integral can be solved using formula 
3.147(7) of 15|, thus leading to the following nontopological pulse-like lattice 

sn'(</>, k) = '\ r , (28) 

Zl + [1 — Zl) sn^(w, sj 

with w = a/2 \p\ k^ Zl (1 — Z2) X and the modulus s given by the same expression as above, 
Eq. (j25p . Of course, one notices immediately that this soliton lattice is different from 
the previous one, Eq. ((211). o'^^ period 2L = 2 K (s) / ^y 2 \p\ k'^ zi {1 — Z2), oscillates 
between sn"-'^^^ and K{k) [respectively, —sn^^y^ and —K{k)]. 

The energy corresponding to one period 2 L of this pulse lattice is given by: 



SL 



2\p\ 



Zl (1 - Z2) 



(1 - Vk'f _ a 



K{s) + n(-^ ^, s 

-21 



A plot of the solutions is given in Fig. 4(b) where the solid and dashed lines correspond, 
respectively, to the Solution 1 and Solution 2. Therefore, as expected, we obtained two 
different soliton lattices (with two different limiting cases of topological single solitons). 
Note that the terms "large" and "small" applied to the topological solitons do not refer, as 
usual, to the length of the interval they cover in space, but to the height of the barrier 
they "encounter" P]. 

Case II.2(ii) : When V}^^^ ^ = \p\e a'^ < Vmax, i.e., (E - k'f/e ^ < 
(1 - Tk'f/k^ there is no soliton lattice solution. 



C. Case III : p > and g 

Unlike the last two cases, here p"^ can be > or < g^. However, it turns out that in this 
case, whatever the value of F = \/\q\/p, the potential Val{4>) has only two extrema in one 
period ^ < 2 K{k), namely an absolute minimum Vmin = for = and a maximum 
Vmax = pk'^{l + F^) for = K{k). Note that we chose C = pfc^F^ in Eq. (P). A plot of 
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K(k) 



K(k) 2K(k) 2 




FIG. 4: Case II. 2 : k' < T < 1. (a) The shape of the potential Val(</')- The soHd horizontal line 
corresponds to the integration parameter A^. (b) The two nontopological soliton lattice solutions 
[i.e., the field (/>(x)] corresponding to the value of in (a); the solid line represents Solution 1 (see 
the main text), while the dashed line-Solution 2. 



the potential Val(0) is given in Fig. 5(a). With z = sn^(0, k), Eq. © now becomes: 



2^/2pk^x = ± / 



dz 



with 



2^1,2 



'z y/{zi - z){l- z) {z - Z2) {z - 0) 
;i + k'^ + e) 



2k^ 



1 ± Wl 



:i + F2A;'2 + a^k^y 



(30) 



(31) 
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zi>l^z>Z2>0. Using the formula 3.147(5) of Q|, one can evaluate the integral in 
Eq. thus obtaining a nonto;^o/o5'ica/ soliton lattice 

[(^1 - Z2) - 2i (1 - Z2) sn2(?/, t)] 



sn (0, k) 



[{z, - Z2) - (1 - Z2)^n\y,t)] 



(32) 



Here y = ^J2pk'^ {zi — Z2) x and the elliptic modulus: 



< t 



'^ 1(1 - Z2) 

Zl - Z2 



^ 1. 



(33) 



In one period 2L = 2 K{t) / ^J2pk'^ {zi — Z2), (f) oscillates between sn^^^/z2 and K{k) 
[respectively, —K{k) and — sn~^y^]. A plot of the solution is given in Fig. 5(b). 
The energy corresponding to one period 2 L of this lattice is: 



EsL 



2p 



Zl - Z2 

Zl - 1 



Zl Kit) - {zi - 1) n 

(1 - Z2){Zik^ 







\Zl - Z2 ) _ 





Zl 



zik^ - 1 



m- 



(zi ^2 



n 



A;'2 {zi 



Z2) 



(34) 



IV. SINGLE SOLITON SOLUTIONS OF THE AL POTENTIAL 



Having obtained the soliton lattice solutions of the AL problem, it is now straightforward 
to consider the specific limits of the integration constant and to obtain the corresponding 
single soliton solutions of the Lame potential and compute their energy. 

Cases I.l and I.2(i) : The limit of a single soliton (that corresponds to the period of 
the lattice L 00) is obtained for t / 1 [a^ \ 0; ;zi ^ (1 - T^k'^)/^ and Z2 \ 0]. 
One notices that now (p can vary continuously between —K{k) and K{k) while x explores 
the whole real axis; therefore, one obtains a topological single- soliton (kink-like) solution. 



sn(0, k) 



T^k'^ 



(35) 



1 - k'^ + k'^{l - r^)sinh\y* 

where y* = \/2pk'^ {\ — k"^) x. The corresponding energy of this single kink-like soliton 
is obtained from Eq. (fT^ . using the relations 



lim^^iH(n^r) = lim^^i- - — 

1 — 2(1 — 77,2) 



In 



1 — n 



(36) 
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K(k) 



K(k) 2K(k) 




FIG. 5: Case III : p > and q ^ 0. (a) The shape of the potential VAiifp)- The sohd horizontal 
hne corresponds to the integration parameter A"^. (b) The sohton lattice solution [i.e., the field 
0(a;)] corresponding to the value of in (a). 

It may be noted that all the divergences cancel mutually, thus leading to a finite energy 
of the single soliton: 

'VT 



E 



K 



1, IVl-T^k'^ + k 



PA;'^ - k 



Vi - m'^ + (Tk) 

(r k) \ Vl - T^k'^ - (Tk) 



In 



(37) 



Following Manton 



12| one can compute the asymptotic interaction energy between two 
solitons in the array, or between a soliton and an anti-soliton. As an illustration we will 
consider this particular soliton solution; for all the other cases one can follow exactly the 
same procedure. 
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'4(1 -r^A;'^)" 




exp ( 


_ k'^ (1 - r2) _ 



From Eq. one obtains the asymptotic shape of the sohton, e.g., for x — > +00: 

exp(^-^/2jW{l^^Wk^x^ . (38) 

Then, if 2 L ^ 1 is the distance between two sohtons in the array, according to ^IJj the 
asymptotic interaction energy is 

f/(2 L) ^ ^ exp (-2 L ^2^ P (1 - P fc'^)) . (39) 

If one considers now the small parameter ^ 1 that measures the distance from the single 
soliton limit, according to the results in Sec. Ill we can obtain the asymptotic expression 
for 2 L as: 

1 /16 fl - A;'^)^\ 

and thus 



that corresponds to a repulsive asymptotic interaction between two solitons in the array. One 
can consider also the asymptotic interaction energy between a soliton and an anti-soliton 
and obtains simply U{a^) ~ ~\/ 2{i-r^ k''^) attractive interaction. 



Sine-Gordon Limit 

Taking now the limit k —>■ 0, with \p\, \q\ 00, so that \p\ k"^ ^ P = finite and 
\q\k'^ Q = finite, in the expressions of the AL single soliton solution, Eq. ()35|). one 
obtains the correct kink-like solution of the sine-Gordon potential in Eq. namely 

sin0,G = ±sech(y*) , (42) 



with y* = a/2P(1 — F^) x, and the energy 



EsG = 2v/2P(l - F2). (43) 

Case I.2(ii) : The single soliton limit L ^ 00 [i.e., n y 1; \ (1 — F^), zi — > 1, 
and Z2 {I — r^)/^^] reads: 



1 -F2 

1 - F2 + (F2 - A;'2)sech'(w*) 



sn(0. A;) = ± J- , ,.,^_,„,2, ..^ ' (44) 
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with u* = a/2p (F^ — k"^) x. It represents a nontopological soliton, with (f) varying 
between K{k) and sn^^(A/(l — T'^)/k'^) [respectively, —K{k) and —s'n^^{^J{\ — V'^)/k'^)] 
when X runs over the real axis, i.e., it does not connect two adjacent degenerate minima of 
the potential. 

Its energy is obtained from Eq. (fT^ and reads: 



2pF 



In 



{Tk) + VP - k'^ 



- -In 



k + VP - A;'2 



(45) 



Case II. 1 : Consider now the single soliton limit r 1 [a^ \ 0, Zi 
and 2:2 ^1]- One obtains the following topological (i.e., kink-like) soliton: 



sn(0. A;) = ± 



^'2 _ r2 + (1 _ r2)sinh2(?/*; 

with = y^2~[p["P~(^72"^^^T2y X, of energy 



sinh(?/* 



;i - v^)/k\ 



(46) 



1 / 

— arctan = 
k VVA;'2 - r2 



(FA:) 



arctan 



(rfc) 



(47) 



where we took the limit r 1 in Eq. (PT|) using the relation 

K{r 



Ymvr /iTli—n ,r) = lim,. 



/r 



+ 



arctan(n) . 



(48) 



1 + n2 1 + 71,2 

Note that one can also consider the limit of the sine-Gordon potential in this case and 
obtain the correct sine-Gordon soliton. 

Case II.2(i) : Solution 1 : Large Topological Kink : Considering the single soliton 
limit, s 1 [a2 \ 0, Zi^2 ^ (1 ~ rA;')/A;2], one obtains a "large" topological kink, that 
interpolates between two adjacent minima "across" the large barrier Vmax of the potential: 



sn(0, k) = ± 



1 - Tk' 



k' (r - k') + (1 - rk')tanh\u*) 



tanh(ci;*) , 



(49) 



with uj* = ^2\p\ k'{l - Tk') (E - k') x. Its energy is found from Eqs. ()2fij) and ()48j) as: 



E 



K 



2 V2N 



1 - FA;' ^ /A;'fl - FA;') 

arctan^ / - — — — — F arctan-i 



k' (F 



F - A;' 



(50) 



Case II.2(i) : Solution 2 : Small Topological Kink : Considering now the single 



soliton limit s 1 [a^ ^ 0, 



1 — FA;')/A; ], one obtains a "small" topological 
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kink, that interpolates between two adjacent minima "across" the small barrier V^^^ of the 
potential: 

I 1-Tk' 
' ~ \J 1 - Tk' + k'{T - k')tanh\iu*) ' ^ ^ 

with uj* = ^2\p\ k'{l - Tk') {T - k') x. Its energy is found from Eqs. (j29j) and (j48j) : 



E 



K 



2^2 



r arctan^ 



k' 



arctan 



k' (r - k') 



(52) 



k'{l - Tk') S I - Tk' 

Case III : The single soliton limit t /" 1 [a^ \ 0, zi ^ (1 + r^fc'^)/^^ and Z2 0] 
represents a topological soliton 



sn(0, /c) = ± 



1 + r2 A;'2 + A;'2(l + r2)sinh2(?/*) 
with ?/* = a/2p (1 + A;'^) x. Its energy is given by: 



(53) 



E 



K 



^/2pk^ 



1 / Vl + r2 A;'2 + A; 
A? ^ Wl + r2F2 - A; 



2P 



(FA:) 



arctan 



:rA:) 



(54) 



Vl + r2A;'2 

Here again one can consider the limit of the sine-Gordon potential and obtain the correct 
sine- Gordon soliton of energy Esg = 2 -\/2 P (1 + r2) . 



V. SOLITON LATTICE AND SINGLE SOLITON SOLUTIONS OF THE LAME 
POTENTIAL 



We shall now demonstrate that having obtained the AL soliton lattice and single soliton 
solutions, one can immediately obtain the corresponding solutions of the Lame problem by 
taking suitable limits of the corresponding AL solutions. 



A. Lame Soliton Lattice Solutions 



Starting from the AL soliton lattice solutions and considering the limit F — 0, we obtain 
the two following types of Lame soliton lattice solutions: 

Cases I.l and III {p > 0) : These lead to type I Lame soliton lattice solution. The 
solution is given by Eq. Q with simpler forms for zi,Z2, namely zi = l/k"^, Z2 = a?. As a 
result the Lame soliton lattice energy has the simpler form 

a''K{t) 



EsL 



2p 

Zl - Z2 



z^Kit) - (zi - i)n 



Z2 

Z2' 



(55) 
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Case II. 1 (p < 0) : This produces type II Lame soliton lattice solution. In this case 
the solution is in fact the same as that given by Eq. (fTUj) . except that now zi, take the 
simpler values Z\ = Xjk^ ^ Z2 = ^ — a?- As a result the Lame soliton lattice energy has the 
simpler form 



E, 



SL 



2\p\ 



Z\ - Z2 



-{z, - l)K{r) + ^iH - 



Z2 



Zl - Z2 



a^K{r) 



(56) 



B. Lame Single Soliton Solutions 

We can now easily obtain the corresponding Lame soliton solutions, either by taking the 
appropriate limits of the Lame soliton lattice solutions, or of the AL soliton solutions. 

Type 1 {p > 0) : For example, considering F = in the expressions for the AL soliton 
solution as given by Eq. (jHSl), the corresponding Lame one soliton solution turns out to be 

^ , (57) 



sn((/), k) = —= 

Jl + k'"^ sinh^(?/*) 

where y* = ^2pk^ x. The corresponding energy turns out to be 



E 



2p In 



l + k 
1 - k 



(58) 



where we have used Eq. (jHTj). 

Type II < 0) : Considering F = in the expressions for the AL soliton solution as 
given by Eq. the corresponding Lame soliton solution turns out to be 

sinh(?/*) 



sn((/), k) = ■ 

\ k''^ + sinh^(?/*) 

where y* = \p2\p\WkP' x. The corresponding energy turns out to be 

E = 2^/2\p\ arctan (^^^ , 

where we have used Eq. (jlTj). 



(59) 



(60) 



VI. CONCLUSION 



We have presented the exact soliton lattice and single so. 
corresponding energies, for the Associated Lame equation 



iton solutions, as well as their 



a 121, 110(1 in various parameter 
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regimes. In appropriate limits we also obtained the single soliton and soliton lattice solutions 
of the Lame equation. The topological and nontopolopcal nature of the different solutions 
was discussed. As an illustration of Manton's method |12] we also computed, in a particular 
case, the asymptotic interaction energy between these solitons. In addition to their relevance 
in the study of nonlinear phenomena, these solutions provide valuable information about 
domain walls in field theory, materials and many physical systems. It would be worthwhile 
to analytically study the exact thermodynamical properties of these systems. 
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